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Abstract. The production of prompt photons at high transverse momentum in proton-proton 
collisions (p-p) is a useful tool to study perturbative Quantum-Chromo-Dynamics (pQCD). In 
particular, they yield valuable information about parton distribution functions in the proton. 
The experimental measurement of prompt 7 is a difficult task due to the large background of 
decay photons from neutral mesons, mainly n°. We present a full simulation and reconstruction 
study of prompt 7 identification in p-p at ^/s — 14 TeV in the ALICE electromagnetic 
calorimeter EMCal, giving details on the methods developed to separate them from decay 
photons with the help of shower-shape and isolation cuts. We present Monte Carlo predictions 
for signal and background. The method used to extract the final isolated 7 corrected cross- 
section is presented and the calculation of various experimental corrections is outlined. 



1. Physics motivations 

The study of photon production at large transverse momenta (px 3> A-qcd = 0.2 GeV) in 
hadronic interactions is a valuable testing ground of the perturbative regime of Quantum 
Chromodynamics (pQCD) pQ. At the LHC, photons will allow one to confront the data with 
pQCD predictions at energies never reached before. Since these photons come directly from 
parton-parton hard scatterings, they allow one to constrain the gluon distribution function in 
the proton at small parton momentum fraction x = p pa rton/Pproton- Also, photons produced in 
p-p collisions provide a vacuum baseline reference for the study of their production rates in 
nucleus- nucleus collisions. 

At lowest order, three partonic mechanisms produce prompt 7 in hadronic collisions: (i) quark- 
gluon Compton scattering qg — > 7(7, (ii) quark-antiquark annihilation qq — > 75, and (iii) the 
collinear fragmentation of a final-state parton into a photon, e.g. qq —> qq —> ~fX. The photons 
produced in the two first point-like processes are called direct, the latter fragmentation photons. 
The Compton channel is particularly interesting as it provides information on the proton gluon 
distribution [2] , which is otherwise only indirectly constrained via the "scaling violations" of the 
proton structure function in deep-inelastic scattering (DIS) e-p collisions [3]. 

The measurement of high-p^ photon production is complicated by a large 7 background 
from hadrons, specially from 7r° mesons, which decay into two photons. At high px there are 
between 10 and 100 (at 100 GeV/c and 10 GeV/c respectively) times more ir° than prompt 
7. Furthermore, above 30 GeV/c, the two ir° decay photons merge into a single cluster in 
EMCal and thus cannot be identified as a 7r° via 7-7 invariant mass analysis. In order to 
measure the prompt 7 signal out of the overwhelming background of 7r°, which are produced in 



the fragmentation of jets, one requires the photon candidate to be isolated from any hadronic 
activity within a given distance around its direction. The corresponding measurements are then 
dubbed isolated photons. A standard isolation requirement is that within a cone around the 7 

direction defined in pseudo-rapidity 77 and azimuthal angle (f> by R = \J (v ~ f]- t ) 2 + (</> — (t>/) 2 -, 
the accompanying hadronic transverse energy is less than a fixed fraction e (e.g. often 10%) of 
the photon's pr- R is usually taken between 0.4 and 0.7. Isolation enables one to reject a large 
fraction of the 7r° decay photons, without supressing too many photons of interest. The bigger R 
and the smaller e, the more selective is the isolation. But the choice of optimal isolation criteria 
also needs to consider detector acceptance and Underlying Event (UE) limitations. The values 
choosen are R = 0.2 (for the first year of running at y/s = 7 TeV, with only 4 Super-Modules 
of EMCal installed), R = 0.4 and e = 0.1 (which e.g. removes UE-hadrons with px <2 GeV/c 
in events with a 20-GeV 7). 

2. Monte Carlo predictions 

The JETPH0X [4 J program (for JET-PHOton/hadron X-sections) allows one to calculate 
perturbatively the 7 cross-sections at Next to Leading Order accuracy. It includes both the 
fragmentation and the direct components and isolation cuts can be applied at the parton 
level. Figure [T] shows the expected isolated photons (R =0.2, e=0.1) cross-sections at 7 TeV 
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Figure 1. Differential cross-section as a function of px for isolated photons (R =0.2, e=0.1) in 
p-p collisions at y/s = 7, 10 and 14 TeV (left) in one unit of rapidity at mid-rapidity and expected 
isolated yields and pt reach at y/s = 7 TeV, for an integrated luminosity of L = 15 pb _1 (right). 

and, for comparison, at 10 and 14 TeV (left) and the corresponding yields in p-p at 7 TeV 
(right). The expected pt reach, which will be accessible in the ALICE EMCAL acceptance 
(At/ = 1.4, A0=1.1), during the first year of LHC running at y/s = 7 TeV, considering an 
integrated luminosity of L = 15 pb _1 , is found to be around 80 GeV/c. PYTHIA 6.420 [5J, with 
the Perugia hard and soft tunes [BJ, has been used to study the signal and the background, for 
p-p collisions at y/s = 14 TeV. We have generated the different photon subprocesses: Compton 



We assume fully-efficient high-pr EMCal cluster triggers with no scale-down reduction. 
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Figure 2. Fraction of isolated (R = 0.4, e = 0.1) over inclusive photons and tt° for p-p collisions 
simulated by PYTHIA at y/s =14 TeV (left), and signal (7) over background (tt ) after isolation 
cuts for both the FF hard and soft PYTHIA tunes (right). 



(MSUB = 29), annihilation (MSUB = 14) and fragmentation (MSEL = 14), as well as tt° (MSEL = l) 
events, over one unit of rapidity at y = 0, for different pr bins: [5-20], [20-50], [50-100], [100- 
250], [250-500] and [500-1000] GeV/c, with 10 millions events in each. Isolation cuts (R = 0.4, 
e = 0.1) have been applied on them. Figure [2] (left) shows the fraction of isolated over inclusive 
particles, for the different photon subprocesses and for 7r°, at the MC level. For the direct 7, 
which are theoretically all isolated, we have studied the effect of the UE on the isolation by 
using two different PYTHIA tunes, hard and soft, which differ by their proportions of initial and 
final-state radiation, multiple interactions, beam remnants and parton-to-hadron fragmentation 
functions (FF). The fraction of isolated direct photons (R = 0.4, e = 0.1) goes from 80% to 
100% from 10 to 100 GeV. Thus, the UE leads to a maximum of 20% (±5%) of isolated direct 
photons loss at 10 GeV. The fraction of isolated fragmentation photons with PYTHIA is around 
60%. This leads to a total amount of isolated over inclusive 7 of about 70-80% over the full 
Et range. The isolated tt° fraction is represented by high-z 7r° which carry a large fraction of 
the jet energy (z = Phadron/ Pparton) and are thus isolated from accompanying hadronic activity. 
With PYTHIA we find a fraction of isolated tt° going from 20% to 3% from 10 GeV to 60 GeV. 
High-z isolated 7r° are thus considered as the main background in our prompt isolated photon 
analysis. 

Figure [2] (right) shows the MC-level predictions for the signal over background (S/B) with 
isolation cuts, for both PYTHIA tunes, in order to assess the systematic uncertainty linked to the 
choice of the FF. Isolation allows one to strongly enhance the signal over background ratio. Yet, 
with isolation cuts, the S/B remains smaller than 1 up to ~ 30 GeV/c. In order to suppress the 
remaining background due to isolated 7r°, one thus needs to employ shower-shape cuts which 
improve the photon over ir° ratio, and subtract statistically any final remaining 7r° background. 

3. Analysis 

The goal of this analysis is to obtain a fully corrected isolated prompt photon spectrum at 
14 TeV. The isolated photon cross-section can be written as : 

gd?a 1 1 /N (clusters iso. et id. 7) da n ^o 

dp 3 2-irpT p(t|7) \ A x L x e x Apx x A?? dpx 




with N(clusters id. 7 and iso.) the number of clusters identified as photons and isolated, which 
represents the Signal+Background, p(7|7) the gamma identification efficiency and p(7|vr°) the 
identification contamination due to ir° identified as photon, L the integrated luminosity, A the 
acceptance term, and e all the other corrections (reconstruction and UE isolation, shown in 
Fig. [2] left). For the acceptance term, we consider a fiducial cut corresponding to an isolation 
radius of R = 0.4, which represents 1/4 of EMCal acceptance. The acceptance with respect to full 
coverage for one unit of rapidity at mid-rapidity is thus found to be around 10%. Full simulation 
and reconstruction events for signal (7-jet) and background (di-jets) at 14 TeV have been used 
to obtain all correction factors in Eq. (fT]). The clusters are identified by the shower shape 




Figure 3. PID contamination of the 7r° identified as 7, p(7|7r°) (left), and ratio data/theory 
(full sim reco vs PYTHIA) for isolated 7 (R = 0.4, e = 0.1) in p-p at y/s =14 TeV (right). 



Bayesian method [7]. A track-matching cluster-charged particle cut, using information from the 
central tracker, is also applied on the photon candidates. The photon identification efficiency, 
^(7)7), is found to be constant at ~80% up to 45 GeV/c. The background contamination due 
to 7r°'s identified as 7, p(7|vr°), is found to be less than 10% above 10 GeV/c and decreasing 
withpr (Fig. _J|left). Our shower-shape cuts remove, thus, a significant fraction of neutral pions. 
The reconstruction efficiency is ~90% up to ~80 GeV/c. The signal/background after photon 
reconstruction, identification and isolation is found to be greater than one above pt ~ 15 GeV/c. 

Finally, we use the PYTHIA spectrum of isolated ir° to statistically remove any remaining 
contamination. This 7r° cross-section, weighted by the term p(7|vr°), is subtracted from the 
fully-corrected left term of Eq. (1) to obtain the final corrected isolated 7 spectrum at 14 TeV. 
The systematic errors propagated on the final yields are dominated by the EMCal energy scale 
(about 15%) and the (theoretical) uncertainty in the subtracted isolated- 7r°'s (~20%). Fig. [3] 
(right) shows the final data/theory ratio, i.e. the ratio of the full sim+reco spectrum over the 
original PYTHIA distribution. The fully corrected spectrum is in good agreement with the input 
MC LO spectrum within the expected statistical (for 24 pb -1 ) and systematic uncertainties. 
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2 Defined as a conditionnal probability of identifying a reconstructed cluster as a 7 knowing that it is indeed a 7. 



